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ABSTRACT 

Type la supernovae (SNe la) play an important role in the study of cosmic evolution, 
especially in cosmology. There are several progenitor models for SNe la proposed in 
the past years. In this paper, by considering the effect of accretion disk instability on 
the evolution of white dwarf (WD) binaries, we performed detailed binary evolution 
calculations for the WD + red-giant (RG) channel of SNe la, in which a carbon-oxygen 
WD accretes material from a RG star to increase its mass to the Chandrasekhar mass 
limit. According to these calculations, we mapped out the initial and final parameters 
for SNe la in the orbital period-secondary mass (logP 1 — M%) plane for various WD 
masses for this channel. We discussed the influence of the variation of the duty cycle 
value on the regions for producing SNe la. Similar to previous studies, this work also 
indicates that the long-period dwarf novae offer a possible ways for producing SNe 
la. Meanwhile, we find that the surviving companion stars from this channel have a 
low mass after SN explosion, which may provide a means for the formation of the 
population of single low-mass WDs (<O.45M ). 
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1 INTRODUCTION 

Type la Supernovae (SNe la) are excellent cosmological dis- 
tance indicators due to their high luminosities and remark- 
IL* ■ able uniformity. They have been applied successfully in de- 
. £^ | termining cosmological parameters (e.g. f2 and A; Riess et 
■ al. 1998; Perlmutter et al. 1999). However, several key issues 
•_h ' related to the nature of their progenitors and the physics 
of the explosion mechanisms are still not well understood 
(Hillebrandt & Niemeyer 2000; Podsiadlowski et al. 2008; 
Wang et al. 2008a), and no SN la progenitor system has 
been conclusively identified from before the explosion. It 
is generally believed that SNe la are thermonuclear explo- 
sions of carbon-oxygen white dwarfs (CO WDs) in bina- 
ries (for the review see Nomoto et al. 1997). Over the past 
few decades, two families of SN la progenitor models have 
been proposed, i.e. the double-degenerate (DD) and single- 
degenerate (SD) models. Of these two models, the SD model 
is widely accepted at present. It is suggested that the DD 
model, which involves the merger of two CO WDs (Iben & 
Tutukov 1984; Webbink 1984; Han 1998), likely leads to an 
accretion-induced collapse rather than to an SN la (Nomoto 
& Iben 1985). For the SD model, the companion is proba- 
bly a main-sequence (MS) star or a slightly evolved subgiant 
star (WD + MS channel), a red-giant star (WD + RG chan- 
nel), or an He star (WD + He star channel) (e.g. Hachisu et 
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al. 1996; Li & van den Heuvel 1997; Langer et al. 2000; Han 
& Podsiadlowski 2004, 2006; Chen & Li 2007, 2009; Meng 
et al. 2009; Lu et al. 2009; Wang et al. 2009a,b; Wang, Li & 
Han 2010). Note that, some recent observations have indi- 
rectly suggested that at least some SNe la can be produced 
by a variety of different progenitor systems (e.g. Patat et al. 
2007; Voss & Nelemans 2008; Wang et al. 2008b; Justham 
et al. 2009). 

An explosion following the merger of two WDs would 
leave no remnant, while the companion star in the SD model 
would survive and potentially be identifiable. Tycho's super- 
nova is a Galactic SN la. Ruiz-Lapuente et al. (2004) find in 
the remnant region that Tycho G, a star similar to the sun 
but with a lower gravity, moves at more than three times 
the mean velocity of the stars there. They argued that Ty- 
cho G could be the surviving companion of the supernova. 
Note that there has been no conclusive proof yet that any 
individual object is the surviving companion star of a SN la. 
It will be a promising method to test SN la progenitor mod- 
els by identifying their surviving companions (e.g. Wang & 
Han 2009). 

The WD + RG channel is a possible ways to produce 
SNe la, and is supported by some observations. It is sug- 
gested that, RS Oph and T CrB, both recurrent novae are 
probable SN la progenitors and belong to the WD + RG 
channel (e.g. Belczyhski & Mikolajewska 1998; Hachisu et 
al. 1999a, 2007; Sokoloski et al. 2006). Meanwhile, by de- 
tecting Na I absorption lines with low expansion velocities, 
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Patat et al. (2007) suggested that the companion of the pro- 
genitor of SN 2006X may be an early RG star. Additionally, 
Voss & Nelemans (2008) studied the pre-explosion archival 
X-ray images at the position of the recent SN 2007on, and 
they consider its progenitor as a WD + RG system. 

Xu & Li (2009) recently emphasized that the mass- 
transfer through the Roche lobe overflow (RLOF) in the 
evolution of WD binaries may become unstable (at least 
during part of the mass-transfer lifetime). This important 
feature has been ignored in nearly all of the previous the- 
oretical works on SN la progenitors except for King et al. 
(2003) 1 and Xu & Li (2009), who inferred that the mass- 
accretion rate onto the WD during dwarf nova outbursts can 
be sufficiently high to allow steady nuclear burning of the 
accreted matter and growth of the WD mass. Following the 
work of Xu & Li (2009), Wang, Li & Han (2010) studied the 
WD binaries towards SNe la systematically. However, the 
interest of Wang, Li & Han (2010) mainly focused on the 
WD + MS channel of SNe la. 

The purpose of this paper is to study the WD + RG 
channel towards SNe la in a comprehensive manner, and 
to show the final parameter spaces of companions at the 
moment of SN la explosion. In Section 2, we simply describe 
the numerical code for the binary evolution calculations. The 
binary evolutionary results are shown in Section 3. Finally, 
a discussion is given in Section 4. 



2 BINARY EVOLUTION CALCULATIONS 

In our WD binary models, the lobe-filling star is a RG star. 
The star transfers some of its material onto the surface of the 
WD, which increases the mass of the WD as a consequence. 
If the WD grows to 1.378 Mq, we assume that it explodes as 
an SN la. We use Eggleton's stellar evolution code (Eggle- 
ton 1971, 1972, 1973) to calculate the WD binary evolution. 
The code has been updated with the latest input physics 
over the past three decades (Han et al. 1994; Pols et al. 
1995, 1998). RLOF is treated within the code described by 
Han et al. (2000). We set the ratio of mixing length to local 
pressure scale height, a — l/H p , to be 2.0 and set the con- 
vective overshooting parameter, 8ov, to be 0.12 (Pols et al. 
1997; Schroder et al. 1997), which roughly corresponds to 
an overshooting length of ~0.25 pressure scaleheights (Hp). 
The opacity tables are compiled by Chen & Tout (2007). In 
our calculations, we use a typical Pop I composition with H 
abundance X = 0.70, He abundance Y — 0.28 and metallic- 
ity Z = 0.02. 

During the mass-transfer through the RLOF, the ac- 
creting material can form an accretion disk surrounding the 
WD, which may become thermally unstable when the ef- 
fective temperature in the disk falls below the H ionization 
temperature ~6500 K (e.g. van Paradijs 1996; King et al. 
1997; Lasota 2001). This also corresponds to a critical mass- 
transfer rate below which the disk is unstable. Similar to the 
work of Xu & Li (2009), we also set the critical mass-transfer 
rate for a stable accretion disk to be 
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M cr , disk ~ 4.3 x 10~ 9 (P orb /4hr) 1 - 7 M yr- 1 , (1) 

for WD accretors (van Paradijs 1996), where P or b is the or- 
bital period. The locations of various types of cataclysmic 
variable stars (e.g. the UX UMa, U Gem, SU UMa, and 
Z Cam systems) in a (P or b, M crj disk) diagram are well de- 
scribed by this expression (Smak 1983; Osaki 1996; van 
Paradijs 1996). If the mass-transfer rate, I-M2I, is higher 
than the critical value M cr ,disk, we assume that the accre- 
tion disk is stable and the WD accretes smoothly at a rate 
Mace = I-M2I; otherwise the WD accretes only during out- 
bursts and the mass-accretion rate is M acc = \Mz\/d, where 
d is the duty cycle. The mass-accretion rate is M acc = 
during quiescence. King et al. (2003) showed that for typi- 
cal values of the duty cycle ~0.1 to a few 10 -3 the accretion 
rates onto the WD during dwarf nova outbursts can be suffi- 
ciently high to allow steady nuclear burning of the accreted 
matter. The limits on the duty cycles of dwarf nova out- 
bursts come from observations (Warner 1995): (1) The out- 
burst intervals for each object are quasi-periodic, but within 
the dwarf nova family, the intervals can range from days to 
decades. (2) The lifetime of an outburst is typically from 2 
to 20 days and is correlated with the outburst interval. The 
quasi-periodicity of the dwarf nova outbursts allows to use 
a single duty cycle to roughly describe the change in the 
mass-transfer rate, though we note that this is a simplifica- 
tion of the real, complicated processes. Similar to previous 
studies (e.g. Xu & Li 2009), we also set the duty cycle to be 
0.01. Meanwhile, we also do some tests for a higher or lower 
value of the duty cycle in our calculations. 

Instead of solving stellar structure equations of a WD, 
we adopt the prescription of Hachisu et al. (1999b) for the 
mass-growth of a CO WD by accretion of H-rich material 
from its companion. The prescription is given below. If the 
mass-accretion rate of the WD, M acc , is above a critical 
rate, M cr ,WD, we assume that the accreted H steadily burns 
on the surface of the WD and that the H-rich material is 
converted into He at a rate M cr> \yD. The unprocessed matter 
is assumed to be lost from the system as an optically thick 
wind at a mass-loss rate M w i n d = M2I — M cr ,WD- The critical 
mass-accretion rate is 

Mcr.wD = 5.3xl0~ 7 (L7 ~ X) (Mwd/M s -0.4) M yr' 1 ,(2) 
X 

where X is the H mass fraction and Mwd is the mass of 
the accreting WD. The following assumptions are adopted 
when |M aC c| is smaller than M C r,WD- (1) When |M aC c| is less 
than M CF| wd but higher than iM cr ,WD, the H-shell burning 
is steady and no mass is lost from the system. (2) When 
|-Macc| is lower than |M C r,WD but higher than |M cr ,wD, a 
very weak H-shell flash is triggered but no mass is lost from 
the system. (3) When |M acc j is lower than |M ci .,wd, the H- 
shell flash is so strong that no material is accumulated onto 
the surface of the WD. 

We define the mass-growth rate of the He layer under 
the H-shell burning as 

M He = 7?H|M acc |, (3) 

where r/u is the mass-accumulation efficiency for H-shell 
burning. According to the assumptions above, the values 
of r?H are: 
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{Mcr,WD/|Macc|, |Macc| > Mcr.WD, 
1, Mcr.WD > | M acc | > §M cr ,WD, (4) 

0, |Macc| < |Mcr ,WD- 

When the mass of the He layer reaches a certain value, He is 
assumed to be ignited. If He-shell flashes occur, a part of the 
envelope mass is assumed to be blown off. The mass-growth 
rate of WDs in this case is linearly interpolated from a grid 
computed by Kate & Hachisu (2004), where a wide range 
of WD mass and mass-accretion rate were calculated in the 
He-shell flashes. We define the mass-growth rate of the CO 
WD, M C o, as 

M C o = J?h c Mhc = r?Hor?H|Macc|, (5) 

where r/Hc is the mass-accumulation efficiency for He-shell 
flashes. 

The evolution of these WD binaries is driven by the 
nuclear evolution of the donor stars, and the change of the 
orbital angular momentum of the binaries is mainly caused 
by wind mass- loss from the WD. We assume that the mass 
lost from these binaries carries away the same specific orbital 
angular momentum of the WD (the mass-loss in the donor's 
wind is supposed to be negligible, but its effect on the change 
of the orbital angular momentum, i.e. magnetic braking, 
is included). We incorporate the prescriptions above into 
Eggleton's stellar evolution code and follow the evolution of 
the WD + RG systems. We have calculated a large, dense 
model grid, in which the lobe-filling star is a RG star. 



3 BINARY EVOLUTION RESULTS 

3.1 An example of binary evolution calculations 

In Fig. 1, we present an example of binary evolution calcula- 
tions. Panel (a) shows the M2, Mco and Mwd varying with 
time, while panel (b) is the evolutionary track of the donor 
star in the Hertzsprung-Russcll diagram, where the evolu- 
tion of the orbital period is also shown. The binary shown 
in this case is (M 2 , M^ D , log(PVday)) = (1.2, 1.0, 0.2), 
where M\, M WD and P' are the initial mass of the donor 
star and the CO WD in solar masses, and the initial orbital 
period in days, respectively. The donor star fills its Roche 
lobe on the RG stage which results in case B mass-transfer. 
During the whole evolution, the mass-transfer rate is lower 
than the critical value given by Equation (1). Thus, the ac- 
cretion disk experiences instability. The mass-accretion rate 
Macc of the WD exceeds ^M cr ,wD after the onset of RLOF, 
leading to the mass-growth of the WD during outbursts. Af- 
ter about 2.6 x 10 8 yr, the WD grows to 1.378 Mq, which 
explodes as an SN la. At the SN explosion moment, the 
mass of the donor star is M| N = 0.4418 Mq and the orbital 
period log(P SN /day) = 0.7992. 

3.2 Initial and final parameters for SNe la 

In Fig. 2, we show the initial contours for SNe la and the 
final state of binary evolution in the (logP 1 , M 2 ) plane at 
the moment of SN la explosion with the duty cycle d — 0.01. 
The initial WD masses are 1.0, 1.1 and 1.2 M©. The final 
contour (solid line in Fig. 2) is much lower than that of 
the initial contours, which results from the mass-transfer 
from the secondary to the WD and the mass-loss from the 
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Figure 1. An example of binary evolution calculations. In panel 
(a), the solid and dash-dotted curves show M2 and Mwd varying 
with time, respectively. The open circles represent the evolution 
of Mco during outbursts. In panel (b), the evolutionary track 
of the donor star is shown as a solid curve and the evolution of 
orbital period is shown as a dash-dotted curve. Dotted vertical 
lines in both panels and asterisks in panel (b) indicate the position 
where the WD is expected to explode as a SN la. The initial 
binary parameters and the parameters at the moment of the SN 
la explosion arc also given in these two panels. 
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Figure 2. Parameter regions producing SNe la in the orbital 
period— secondary mass plane (logP 1 , Mj) for the WD + RG 
systems with the duty cycle d = 0.01. The initial WD masses are 
1.0, 1.1 and 1.2 M Q . The final state of the WD + RG systems in 
the plane is encircled by the solid line, where filled stars, circles 
and triangles denote that the WD explodes as an SN la with 
initial WD masses of 1.0, 1.1 and 1.2 Mq, respectively. 
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Figure 3. Regions in the initial orbital period-secondary mass 
plane (log P 1 , ) for the WD + RG channel that produce SNe 
la with initial WD mass of 1.2 Mq, but for different duty cycle 
values. 

system. Note that, the enclosed region almost vanishes for 
Mjyu = 1.0 Mq, which is then assumed to be the minimum 
WD mass for producing SNe la from this channel. If the 
initial parameters of a WD + RG system are located in the 
initial contours, an SN la is then assumed to be produced. 
These initial contours for the WD + RG channel can be 
expediently used in the binary population synthesis studies. 
The data points and the interpolation FORTRAN code for 
these contours can be supplied on request by contacting BW. 



4 DISCUSSION 

Compared with most of previous investigations on the WD 
+ RG channel of SNe la, the main difference is that we 
considered the effect of the accretion disk instability on the 
evolution of WD binaries. In our results, there is no WD + 
RG system with period as long as ~10 3 days as indicated 
in Li & van den Heuvel (1997) and Hachisu et al. (1999a). 
This is because if the initial period of a WD binary is too 
long, the mass-transfer rate between the RG donor star and 
the WD will be too high and the optically thick wind will 
occur and take much H-rich material away from the binary 
system. Finally, the RG donor star has no enough material 
to be accumulated onto the WD. 

Similar to previous studies (e.g. King et al. 2003; Xu 
& Li 2009), our study also indicates that the long-period 
dwarf novae offer a possible ways for producing SNe la. In 
particular, there is an advantage for this work, i.e. the SN 
la explosion is always occur in a WD binary of small sec- 
ondary/primary mass ratio (<1), and with very little of the 
H envelope of the secondary remaining in our investigations. 
This feature will greatly reduce the possibility of H contam- 
ination of the SN la ejecta (see also King et al. 2003), which 
is consistent with the defining characteristic of most SNe la 
having no detectable H (Branch et al. 1995). This also has 
broad implications towards the work that is currently plac- 
ing very tight limits on the amount of H entrained in the 
SN la ejecta, determined through the analysis of late-time 
spectra (e.g. Mattila et al. 2005; Leonard 2007). 



However, the results in this paper depend on many un- 
certain input parameters, in particular for the duty cycle 
which is poorly known. The main uncertainties lie in the 
facts that it varies from one binary system to another and 
may evolve with the orbital periods and mass-transfer rates 
(e.g. Lasota 2001; Xu & Li 2009). This is the reason why 
we choose an intermediate value (0.01) rather than other 
extreme values (e.g. 0.1 or 10 -3 ). Furthermore, we also did 
some tests for a higher or lower value of the duty cycle. In 
Fig. 3, we show the influence of the variation of the duty 
cycle value on the regions for producing SNe la with initial 
WD mass of 1.2 Mq. We see that, for a high value (0.1), 
the right boundaries of the regions will be shifted to higher 
period, while, for a low value (10 ), the right boundaries of 
the regions will be shifted to lower period. For the low value 
of the duty cycle, it will have a high mass-accretion rate of 
WDs during outbursts, so that the accreting WDs will lose 
too much mass via the optically thick wind, preventing them 
increasing masses to the Ch mass. Thus, a low value of the 
duty cycle will reduce the regions for producing SNe la. 

In this paper, we set the metallicity Z — 0.02. For 
the WD + RG channel, varying the metallicity would have 
strong influence on the regions for producing SNe la (Fig. 
2), e.g. high metallicity leads to larger radii of zero-age MS 
(ZAMS) stars, then the left boundaries of the regions will be 
shifted to longer period. Meanwhile, stars with high metal- 
licity evolve in a way similar to those with low metallicity 
but less mass (Umeda et al. 1999; Chen & Tout 2007; Meng 
et al. 2009). Thus, for the WD binary systems with par- 
ticular orbital periods, the companion mass increases with 
metallicity. 

At present, the existence of a population of single low- 
mass (<0.45M Q ) WDs (LMWDs) is supported by some ob- 
servations (e.g. Maxted et al. 2000; Kilic et al. 2007). The 
formation of single LMWDs is still unclear. It is suggested 
that single LMWDs could result from single old metal-rich 
stars which experiences severe mass loss prior to the He flash 
(Kalirai et al. 2007; Kilic et al. 2007). However, the study 
of initial-final mass relation for stars by Han et al. (1994) 
implies that only LMWDs with masses > 0.4 Mq might be 
produced through such a single-star channel, even at high 
metallicity (Meng, Chen & Han 2008). Thus, it would be dif- 
ficult to conclude that single stars can produce the LMWDs. 

The companion in the WD + RG channel would survive 
and evolve to a WD finally. In Fig. 2, we see that the com- 
panion stars have a low mass at the moment of SN explosion. 
The companion stars will be stripped of some mass due to 
the impact of SN ejecta. Marietta et al. (2000) presented sev- 
eral high-resolution two-dimensional numerical simulations 
of the impact of SN la explosion with companions. They 
find that a RG donor star will lose almost its entire enve- 
lope (96%— 98%) owing to the impact of the SN la explosion 
and leave only the core of the star. Thus, the surviving com- 
panion stars from this channel will have a relatively low mass 
after SN explosion and evolve to a WD finally, which pro- 
vides a possible pathway for the formation of the population 
of single LMWDs (<O.45M ). Meanwhile, we also suggest 
that the observed, apparently single LMWDs may provide 
evidence that at least some SN la explosions have occurred 
with non-degenerate donor stars (such as RG donor stars). 

The CO WD + RG systems can be formed by binary 
evolution. Wang, Li & Han (2010) find that there is one 
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channel which can form CO WD + RG systems and then 
produce SNe la. In the detailed binary evolution procedure, 
the primordial primary first fills its Roche lobe at the ther- 
mal pulsing asymptotic giant branch stage. A common en- 
velope is then easily formed owing to dynamically unstable 
mass-transfer during the RLOF stage. After the common en- 
velope ejection, the primordial primary becomes a CO WD, 
then a CO WD + MS system is produced. The MS compan- 
ion star continues to evolve until the RG stage, i.e. a CO 
WD + RG system is formed. For the CO WD + RG systems, 
SN la explosions occur for the ranges Mi,i ~ 5.0— 6.5 Mq, 
M 2ji ~ 1.0-1.5 M Q , and P i > 1500 days, where M u , M 2 ,i 
and P 1 are the initial mass of the primary and the secondary 
at ZAMS, and the initial orbital period of a binary system. 

The WD + RG channel has a long delay time from 
the star formation to SN explosion due to the RG donor 
star with low initial masses (<1.5Mq). Thus, this chan- 
nel can contribute to the old population of SNe la implied 
by recent observations (Mannucci et al. 2006; Totani et al. 
2008; Schawinski 2009). The old population of SNe la may 
have an effect on models of galactic chemical evolution, since 
they would return large amounts of iron to the interstellar 
medium much later than previously thought. It may also 
have an impact on cosmology, as they are used as cosmolog- 
ical distance indicators. 
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